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Abstract 

We study the reversal of magnetization in an isotropic ferromagnetic film free 
from charges by exposing it to a circularly polarized electromagnetic (EM) 

^ ' field. The magnetization excitations obtained in the form of line and lump 

00 

lO ' solitons of the completely integrable modified KP-II equation which is derived 

o 

L^ ' using a reductive perturbation method from the set of coupled Landau-Lifshitz 

and Maxwell equations. It is observed that when the polarization of the EM- 

field is reversed followed by a rotation, for every ^-degrees, the magnetization 

g. is reveled. 
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Studies on developing high density magneto-optic data storage and development of ul- 
trafast magnetic recording have attained momentum due to their high technological impli- 
cations in the very recent times. It has been identified both experimentally and theoretically 
that magnetization reversal in magnetic films is one of the very fundamental issues in mag- 
netic data storage [Ij-Q. In conventional magnetic recording the reversing field is applied 
anti-parallel to the direction of magnetization of the medium in which case the reversal 
takes place at the nano-second level. However, much shorter reversal times can be achieved 
through precessional reversal of the magnetization ||l| . The torque developed between the 
magnetic moment of the medium and the external magnetic field will make the magnetic 
moment to precess at the pico-second scale. The precessional motion is governed by the 
Landau-Lifshitz (LL) equation and recent experiments have shown the validity of the equa- 
tion at the pico-second level ^. Thus the switching process or magnetization reversal can 
be understood by solving the Landau-Lifshitz equation of motion for the applied magnetic 
field. In this context static and time dependent (pulse) applied field switching phenomenon 
have been studied in the past [§-[§]. Magnetization reversal studies in the case of a more 
general applied magnetic field that varies both spatially and temporally is also equally in- 
teresting and it needs attention. In the context of magneto-optics, it is very relevant and 
important to consider the applied field as the magnetic field component of the EM-field 
and thus the problem can be formulated in terms of the LL equation coupled with the 
Maxwell equations. The one-dimensional version of this kind of problems have been studied 
recently in the case of isotropic and anisotropic ferromagnets and soliton modes were found 
to represent the magnetization excitations and the EM-field has also been modulated in the 
form of solitons 0-|^. The purpose of the present paper is to examine the reversal process of 
magnetization in a ferromagnetic film by generating soliton modes during precession of mag- 
netization in the presence of an EM-field. By assuming that there are no free electric charges 
in the medium, we solve analytically the LL equation coupled with Maxwell equations in 
two spatial dimensions using a reductive perturbation method. 

The dynamics of magnetization density in an isotropic charge-free ferromagnetic film 



under the influence of an external EM-field can be expressed in terms of the LL equation |T^ 
dtM = M A [V^M + AH], M^ = 1, where M(a;, y, t) = (M^, M^, M^) is the magnetization 
density and H(a;, y, t) = {H^, H^, H^) is the magnetic field component of the electromagnetic 
field. The first term in the right hand side of the LL equation represents the contribution due 
to spin-spin exchange interaction between the nearest neighbours and the term proportional 
to A (= g^B'i g = gyromagnetic ratio, ^b = Bohr magneton) corresponds to the interaction 
between the magnetization of the medium and the magnetic field component of the EM-field 
(Zeeman term). When H is a constant or a time dependent field set along a specific direction, 
it can be transformed away using the transformation M^ = {M^iiM^) = M^e^*-^-°° 
and the dynamics remains the same as in the case without any external field. However, when 
the field H varies spatially, the effect due to the field cannot be transformed away in this 
fashion. It may be noted that in the LL equation we have not included the phenomenological 
Gilbert damping term because we have assumed that the medium does not contain any free 
charges which on interaction with magnetic electrons will introduce damping. 

The interaction between the magnetic field component of the EM-field and matter or 
material medium can be expressed in terms of Maxwell equations [IT]]. In the absence 



of static and moving charges. Maxwell equations can be written as V^H — V(V.H) = -^ 
df [H + M], where V(= xdx + ydy) and V^(= d1 + dy) are the two-dimensional gradient and 
Laplacian operators respectively, and c = ^== is the velocity of propagation of the EMW 
and eo and /iq are the dielectric constant and permeability of the medium respectively. An 
untreated ferromagnetic material has the constitutive relation H = — — M where B(x, y, t) 
= {B^, B^ ,B^) is the magnetic induction. Now using this constitutive relation, the LL and 
Maxwell equations can be rewritten as 



^iM = M A 



A 

V^M + — B 

Ato 



M' = 1, (la) 

V(V ■ M) - V"M . (lb) 



where the operator A = {df — c^V^). Thus, the set of coupled equations ( [Ta| ) and (jTBp 



completely describe the interaction of the EM-field with the ferromagnetic film and the 



dynamics of the associated fields namely magnetization, magnetic induction and magnetic 
field from the EM-field. 

The nonlinear character of the LL equation ( p!aD makes the analysis difficult in the present 
form and hence we try to solve Eqs. ( p!aD and (^B|) using a reductive perturbation method [jl2 . 



Assuming that the plane wave of the EM-field travels along x-direction in the ferromagnetic 
film (xy-plane) we introduce the wave variable ^ = {x — vt) where v is the group velocity. 
We also introduce slow space and time variables through the stretching C = £^^, C = ^"^V ^^'^ 
T = s^t, where e is the small perturbation parameter. The solutions of Eqs.(l) are then 
expanded asymptotically about uniform values 

F(e, C, r) = Fo(e, C, r) + £Fi(e, C, r) + e'F,{^, C, r) + ..., (2) 

where F stands for the magnetic induction B and the magnetization M. We substitute the 
slow variables and the expansions in the component forms of Eqs.(p!a]) and {^^), collect terms 
corresponding to similar powers of e and solve the equations at different orders. At 0(e°) 
from the component forms of Maxwell equations (Eq.(p]B[)), we obtain the results Bq = 0, 
Bq = kMQ, where a = y,z and k = [eo(c^ — v^)]^^ and from the LL equation (Eq.(|T^)) 
after using the above results obtained we get Mq = 0. The results at 0{e^) thus show that 
a planar anisotropy is developed in the magnetization (Mq — Mq), the magnetic induction 
{Bq — Bq) and in the magnetic field {Hq — Hq) space at the lowest order of expansion. At 
0{e^), from Eq.(|l^) we obtain Bf = 0, 5° = fcMf and from Eq.(||) we find that (9^M^ 
= ^M^M^ and d^M^ = -^M^M^. At 0{e^), after using the results from the lower 
orders, we find from Eq.(|^ that i?f = and 

1 



[B^ - kM^] = 

eok 



(3) 



;/ di'd^B^ + c^f dCf d^'dlB, 

and from Eq. (]Ta|) , we get 

a^Mf = - {M^dlM^ - M^dlM^ - A [M^B^ + M^B^q - M^B^ - M^Bl] } . (4) 

We now try to solve the set of equations (^ and (^ to evaluate the complete set of solutions 
at the lowest nonvanishing order. For this it is advantageous to switch to the polar coordinate 



representation. As the results at 0{e^) show that the magnetic field is restricted to the 
(Hq — Hq) plane, we consider a circularly polarized or rotating magnetic field component of 
the EM- field in this plane by considering Hq = (0, sin6, cosd). This makes us also to choose 
Mo = {0,sin9,cos9), where 9 = 9{C^,(,t) is the angle made between the direction of the 
applied EM-field and the uniform magnetization of the film. Thus Eq. (^) in the above polar 
co-ordinate system after using Eq.(||) becomes 



fidp = 37 



+ 



2vA 



eo 



sin9 f di' f di'dlcos9 - cos9 f d^' f d^'dlsin9 
J— 00 J— 00 J— 00 J —00 

/•? [^ , 

sin9 I d^ drCos9 — cos9 / d^ drsin9 
J— 00 J— 00 



(5) 



where /z = (^ — 1) and 7 = |^. Differentiating Eq.(|^) twice with respect to ^, we get 



d, 



m 

d,9 



d^9 = -fi{dl 

+ 37^5 



+ 37 



cos9 / d^'d^cos9 + sin9 I d^'d^sin9 
J —00 J— 00 



sin9 I d^'d'jcos9 — cos9 I d^'d'^sin9 
J —00 J— 00 



(6) 



where F{9) = dr9 +fid^9.. While writing the above equation, we have rescaled r — > 21m'''"' 
We assume that the fields vary very slowly along y-direction when compared to x-direction 
(i.e.), {d'^9 <^ 9^9). Then we substitute Eq.(^ in the resultant equation obtained after 
differentiating Eq.(P) once with respect to ^ and again after successive integration and 
differentiation we finally obtain d^ [d^f —^fJ-pd^f +fJ,d^f = —3'yd'^f -I-37 9|/ -|- d^fd^f 
+ t^|//_oo '^^'^c/ ' "where / = d^9. When 7 = /i, this is equivalent to the completely 
integrable modified Kadomtsev-Petviashvili (MKP) equation 

1 



drf + iidlf-?>ii 



-f%f-d^w + wdj 



0, 



d^w = dj. 



(7) 



Eq.(^ when fi = i and /i = 1 are respectively known as the MKP-I and MKP-II equa- 
tions. However, in our problem fi cannot be imaginary and hence we have only the MKP- 
II equation for our further analysis. Different types of soliton solutions to the MKP-II 
equation such as line soliton, lump soliton and breather have been found using the d- 
dressing and inverse scattering transform (1ST) methods [|l^. However, as the structure 



of breather solution is very rich we are unable to present them here. The general N-line 



soliton solution is obtained using 9-dressing method |]T3| and the simplest line soliton of the 
MKP-11 equation corresponding to A^ = 1 can be explicitly written in the form f{^, (, r) 
= -2{a, - /30^i/ [ai/3?(e-^ - ^e^)(e-^ - 7ie^)], where 2G = (^ - j-^)^ -(^ - j,)C 
— 4(^ — -ij)r +ln2\ „ ^^ \ and 71 = sgn{ „ ^^ ). Here R is the Kernel and tti and Pi are 
arbitrary real constants used in the 1ST analysis. This solution is nonsingular only if 71 < 
and ai, 4- > 0. The line is regular in ^ and ( and is a constant along a particular direction. 

Lumps are rational solutions which normally decay in all directions in the plane. Unlike 
the MKP-1 case, in the case of MKP-11, the rational solutions of the equation are singular. 
For instance, the simplest (A^ = 1) lump soliton of MKP-II equation corresponding to 
N=l can be written as /(^, C, r) = 2ai/[^ ~ (^ ~^ a ~ ^)^]' where ai is an arbitrary real 
constant. This solution describes the uniform motion of two simple poles of opposite signs 
along a line which are parallel to each other with a distance ai and move with equal velocity 
12ar^ 

Using the above line and lump solitons in the relation / = || , ^ and finally after using in 
the relation connecting Mj^, Mg and Mq we can calculate the components of magnetization at 
the lowest existing order. For example, the line soliton of the x-component of magnetization 



at the lowest existing order (Mf) is found to be [13 



Mr = —, — -'^^^-^^^'^^ _, (8) 

Fig. (la) shows a snapshot of Mf-line soliton at r = 1 unit for v = 0.5, k = 0.25, A = 1.0, 
«! = 5.0, Pi = 2.5 and 71 = —0.5. 
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(a) (b) 

FIG. 1. (a) Line soliton of Magnetization (Mf ) at r = 1 unit for v = 0.5, k = 0.25, A = 1.0, 
ai = 5.0, /3i = 2.5 and 71 = —0.5, (b) Lump soliton of magnetization (Mf ) at r = 1 unit for 
V = 0.5, k = 0.25, A = 1.0 and ai = 2.05. 



The lump 1-soliton solution of the x-component of magnetization (Mf ) is obtained as 

2vai 



Mf 



(9) 



In Fig. (lb), a snap shot of the lump soliton of the a;-component of the magnetization for 
V = 0.5, k = 0.25, A = 1.0 and ai = 2.05 at r = 1 is given. In Fig. (lb) the uniform motion 
of the two simple poles. Pi and P2 separated by a distance of 2.05 units and travelling with 
a velocity of 2.855 units can be observed. Similarly, line and lump soliton solitons for the y 
and z components of magnetization can also be constructed. 

The magnetization states observed in the form of line and lump solitons in the previous 
case correspond to the states when the rotating circularly polarized magnetic field component 
of the EM-field is acting on the ferromagnetic film which is magnetized perpendicular to the 
film at the lowest order. It is interesting to find that when the direction of polarization of 
the magnetic field component of the external EM-field is reversed (i.e. if the initial field is 
left circularly polarized, now in the present case it should be right circularly polarized and 
vice versa) and then rotated continuously, for every ^, n = 0, 1,2, ... degrees of rotation, 
the direction of magnetization in the film is reversed. In other words, the above can be 



achieved by transforming v — > —v and 6^6 + n|. In view of these transformations, we 
now have Mq = (0, cos^, sm^). We use this and the above transformations and also the 
results from Eq.(|^) in Eq.(^). And after repeating the same calculations of the previous case 
as found after Eq.(^) we once again obtain the MKP-II equation (|^. Then we calculate 
the line and lump solitons of the x-component of the magnetization Mf using the relation 
connecting Mq , Mq as done before and the results are found to be Mf = 2v{ai — /9i)^7i 
/[kA{aiPf{e~'^ — ^^e'~^)(e~'^ — 716*^)}] for the line soliton and the lump soliton in the form 

2 

Mf = —2vai/ kA{^ ~ (^ ~^ a ~ ^Y}]- ^^ comparing the line and lump solitons for Mf 
in both the cases, it can be observed that there is a sign change in the expression for solitons 
and hence magnetization in the present case is reversed. 




(a) (b) 

FIG. 2. (a) Reversal of line soliton of the magnetization (Mf ) given in Fig. (la) for the same 
values of the parameters, (b) Reversal of lump soliton of the mangnetization (Mf ) given in Fig. (lb) 
for the same values of the parameters. Both are due to change in direction of polarization (left ^-> 
right) followed by the rotation of the magnetic field by ^ degrees. 

The new configurations of magnetization due to reversal in the case of Mf correspond- 
ing to the line soliton and also lump soliton with reference to the two poles Pi and P2 have 
been demonstrated in Figs. (2a) and (2b) respectively. Thus the reversing of the direction 
of polarization of the magnetic field and a rotation by —degrees reverses the original mag- 
netization states given in Figs.(la,b). As the precessional frequency of the magnetic dipole 
moment in ferromagnets is in the pico-seconds scale, the magnetization reversal can take 

8 



place in a faster rate compared to the conventional reversal process that takes place in the 
nano-seconds scale by the reversal of magnetization of the domains. This has also close cor- 
respondence with the recent experimental results found in ||I[ on in-plane magnetized cobalt 
films in which pulse-like time dependent magnetic field in the plane of the film as short 
as two pico-seconds is able to reverse the magnetization. It was found that when the field 
encompasses right angle to the magnetization of the film the reversal can be triggered by 
even very small fields. Also, our results have very close correspondence with the numerical 
analysis of the magnetization reversal found in refs. PJ^ using four pico-seconds magnetic 
field pulse which has been explained based on the LL equation with Gilbert damping also 
to take into account the relaxation.. 

In this paper we have studied the reversal of magnetization in an isotropic charge-free 
ferromagnetic film when an EM- field is applied on it, by solving the coupled Maxwell equa- 
tions and LL equation in two-spatial dimensions using a reductive perturbation method. 
The results show that at the lowest order of perturbation the ferromagnetic film is mag- 
netized normal to the plane of the film. In the next order of perturbation it is found that 
the magnetization is excited and a coherent magnetization structure in the form of line and 
lump solitons of the MKP-II equation is obtained. Interestingly we found that, when the 
direction of polarization of the applied EM-field is reversed and the field is rotated contin- 
uously, the magnetization gets reversed for every addition of ^-degrees. As this effects the 
magnetization reversal via the precessional motion, it will decrease the reversal time to the 
order of pico-seconds than nano-seconds in conventional magnetization reversal processes. 
This phenomenon has very close correspondence with the recent experimental and numer- 
ical observations of magnetization reversal in cobalt film when ultrashort magnetic pulses 
are applied to it. This interesting phenomenon foresee applications in ultrafast magnetic 
recording in future. 
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